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CRAB  ANGLE  ESTIMATION  WITH  TWO-DIMENSIONAL 
CROSS  CORRELATIONS  AND  CENTROIDS 

1.  INTRODUCTION 


In  Che  Spring  of  1982  the  author  was  assigned  the  task  of  determining 
whether  accurate  crab  angle  estimates  could  be  obtained  by  applying  a 
certain  two-dimensional  cross  correlation  technique  to  the  data  described 
in  Ref.  [1].  The  scope  of  the  effort,  and  the  general  nature  of  the 
results  obtained  are  described  in  Section  2.  Hundreds  of  different  types 
of  crab  angle  estimators  were  examined,  and  they  were  all  found  to  yield 
results  which  were  consistently  bad  for  large  crab  angles  (5°  or  more), 
and  at  best  marginal  for  smaller  angles.  The  probable  cause  for  this 
poor  performance  is  discussed  in  Section  2  and  documented  in  Section  3, 
viz.,  a  dissimilarity  between  the  returns  of  the  two  radars  used  in  the 
system  caused  by  range  differentials,  aircraft  roll,  and  the  crab  angle 
itself.  In  Sections  3-5,  we  discuss  attempts  to  improve  accuracy  by 
numerically  adjusting  the  thresholds  on  one  or  both  of  the  radars,  by 
various  data  weighting  and  data  smoothing  schemes,  and  by  schemes  for 
outlier  removal.  A  representative  sample  of  the  results  are  tabulated  in 
Sections  4  and  5.  Although  the  similarity  between  the  the  two  returns 
was  often  enhanced  by  applying  these  methods,  the  performance  of  the 
resulting  estimators  was  never  satisfactory.  The  explicit  formulas  used 
in  the  calculations  are  given  in  an  appendix.  A  summary  and  conclusions 
are  given  in  Section  6. 

2.  NATURE  OF  THE  RESULTS 


The  crab  angle  (yaw  angle)  measurement  system  discussed  in  this 
report  is  described  in  Ref.  [1].  It  consists  of  the  two  X-band  radars 
recommended  in  Ref.  [1],  which  are  located  on  opposite  sides  of  the 
runway.  For  each  emitted  radar  pulse  recordings  are  taken  of  range  and 
amplitude,  and  the  data  for  each  trial  run  therefore  consists  of  two 
2-dimensional  arrays,  AMP^  (N,R)  and  AMP2  (N,R)  where  AMPj  (N,R)  is  the 
(unsigned)  amplitude  of  the  detection  made  on  the  N'th  pulse  at  the  R’th 
range  by  the  I'th  radar.  The  ranges  R  were  measured  in  increments  of  1 
ft.,  and  the  radar  prf  was  1000  Hz.  The  useful  data  was  found  to  be 
within  the  bounds  0<N£700  and  60<R<160,  and  each  array  therefore  consists 
of  70,000  entries. 

Our  task  was  to  cross  correlate  the  two  arrays,  find  the  lagged 
value  AN  of  N  which  corresponds  to  the  peak  correlation,  and  then  use 
this  value  to  estimate  the  crab  angle  6  according  to  the  formula  (cf. 

Ref.  [1],  p.  10) 

0  *  tan“l  v  At  t  (1) 
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where 

V  =  aircraft  speed, 

D  =  distance  between  the  radars  =  210  ft.,  and 

At  =  AN/1000  =  time  lag  corresponding  to  AN. 

Previous  attempts  had  been  made  to  measure  At  by  means  of  a  certain  one¬ 
dimensional  cross  correlation  technique  (Ref.  [1]);  however,  the  results 
were  unsatisfactory,  and  it  was  hoped  that  the  two-dimensional  cross 
correlation  technique  would  meet  the  accuracy  goal  of  +  1°.  This  hope 
was  not  realized,  as  will  be  seen. 

The  basic  idea  motivating  this  method  is  that  the  aircraft  is 
symmetric  about  its  fore  and  aft  axis  so  that  the  two  data  records  AMP^ 
and  AMP2  should  be  approximately  the  same  except  for  the  time  delay  At 
produced  by  the  crab  angle.  (See  Ref.  [1].)  However,  the  two  data 
records  were  found  to  differ  significantly;  namely,  one  radar  always  had 
more  detections  than  the  other,  the  peak  correlations  between  the  two 
records  were  very  small,  and  the  two  auto-correlation  functions  of  the 
records  AMP^  and  AMP2  (each  correlated  with  a  copy  of  itself  at  lagged 
values  of  N  and  R)  were  different.  This  lack  of  similarity  between  the 
two  data  records  was  apparently  caused  by  range  differentials  as  well  as 
the  crab  angle  itself.  It  also  appears  that  the  aircraft  rolled  as  it 
maneuvered  to  produce  a  crab  under  conditions  of  zero  crosswind.  The 
roll  was  on  the  order  of  5°,  and  could  have  contributed  to  the  lack  of 
symmetry  between  the  two  data  records.  These  results  are  documented  in 
the  next  section. 

Attempts  were  made  to  enhance  the  similarity  between  the  records  AMP^ 
and  AMP2  by  numerically  adjusting  the  thresholds  of  one  or  both  of  the 
radars,  and  by  various  data  smoothing  and  weighting  schemes.  These  modi¬ 
fications  are  described,  and  the  results  documented,  in  the  following 
sections.  The  results  were  again  disappointing. 

The  correlation  of  two  70,000-element  arrays  is  a  formidable  task. 

The  correlation  function  of  two  2-dimensional  arrays  is  itself  a  2-dimen¬ 
sional  array,  and  hundreds  of  values  (each  being  the  sum  of  70,000 
products)  must  be  calculated  for  accurate  peak  location  by  means  of  lobe 
smoothing  and  beam  splitting.  Because  of  its  numerical  complexity,  even 
if  the  cross  correlation  technique  has  been  successful  it  would  have  been 
highly  desirable  to  examine  other  techniques  which  are  easier  to  imple¬ 
ment.  The  centroid  is  a  natural  candidate,  which  besides  being  easy  to 
compute,  also  has  a  certain  attractive  stability  property;  viz.,  it  is 
not  very  sensitive  to  data  elisions.  There  are  various  kinds  of  cen¬ 
troids,  corresponding  to  the  different  possibilities  for  thresholding  and 
data  weighting  mentioned  above.  In  addition,  various  schemes  of  outlier 
removal  were  used  to  enhance  the  similarity  between  the  records  AMPj 
and  AMP2-  The  results  obtained  with  centroids  were  no  better  and  no 
worse  than  those  obtained  by  the  correlation  techniques. 

The  various  methods  for  data  thresholding,  smoothing,  and  weighting, 
etc.,  can  be  combined  in  a  large  number  of  ways  to  produce  different  types 
of  crab  angle  estimators.  In  all,  more  than  300  types  were  examined  and 
the  results  were  never  of  acceptable  performance. 


3.  PRELIMINARY  ANALYSIS  OF  THE  DATA 


Our  initial  examination  of  the  data  was  confined  to  four  trial  data 
runs  (Files  3,  5,  6,  and  7)  for  which  the  crab  angles  were  small,  moder¬ 
ate,  and  large,  the  cross-track  angles  were  negligible,  and  for  which 
there  were  no  problems  in  the  recording  system  which  required  special 
treatment  to  reduce  the  data  to  usable  form.  The  "true"  values  of  crab 
angle  and  cross  track  angle  were  obtained  from  the  optical  data  described 
in  Ref.  [1].  The  sign  convention  is  chosen  so  that  the  crab  angle  is 
positive  for  "left"  crab  (aircraft  pointing  to  Radar  2)  and  negative  for 
"right"  crab  (aircraft  pointing  to  Radar  1).  These  values  are  shown  in 
Table  1.  The  column  labelled  "SUM"  is  the  signed  sum  of  the  crab  angle 
and  cross  track  angle,  which  is  the  apparent  crab  angle  that  would  be 
sensed  by  the  radar  system.  None  of  the  methods  produced  acceptable 
results  for  the  large  crab  angle  in  File  5.  If,  however,  the  results 
were  acceptable  or  even  marginal  for  the  other  "initial"  files  3,  6,  and 
7,  they  were  then  tested  against  three  "additional"  files,  labelled  4,  9 
and  10.  The  data  from  the  additional  three  files  is  always  separated 
from  that  of  the  initial  four  files  by  a  horizontal  line.  As  it  turned 
out,  these  additional  files  only  served  to  confirm  the  previous  results. 

TABLE  1 

Values  of  Crab  Angle  and  Cross  Track 

Crab  Cross 

Angle  Track  SUM 

File  3  -  .25°  0°  -  .23° 

File  5  +5.5  -0.2  +5.3 

File  6  +1.1  0  +1.1 

File  7  -1.4  0  -1.4 


File  4  -2.0  -  .6  -2.6 
File  9  -4.6  -1.4  -6.0 
File  10  +1.8  +  .5  +2.3 


We  now  focus  on  the  lack  of  similarity  between  the  records  AMP^  and 
AMP2.  Table  2  shows  the  total  number  of  detections  ND(1)  and  ND(2) 
made  by  Radars  1  and  2,  and  the  ratio  ND(2)/ND(1).  There  are  six  sets  of 
data,  labelled  A-F.  The  first  set  "A"  is  derived  from  the  original 
untouched  data,  and  shows  that  Radar  2  consistently  made  more  detections 
than  Radar  1.  This  disparity  was  probably  caused  by  the  fact  that  the 
aircraft  was  30-40  ft.  closer  to  Radar  2  than  to  Radar  1.  The  sets 
labelled  B-F  show  the  number  of  detections  that  remain  after  the  data  has 


been  "reduced"  in  various  ways  in  an  attempt  to  enhance  the  similarity 
between  the  two  records.  In  the  sets  B-E  the  outliers  have  been  removed 
by  calculating  the  means  y  and  second  moments  o  of  each  of  the  records 
AMPj(N,R)  in  the  time  (or  pulse  number  N)  direction,  and  then  rejecting 
the  points  (N,R)  whose  N-coordinate  lies  outside  the  two  limits  L  «  p  +  ko, 
where  k  is  some  fixed  parameter.  In  the  results  shown,  k  =  1  or  2.  By  a 
simple  modification  of  a  DO  loop  parameter,  the  reduced  data  set  can  then 
be  re-cycled  through  the  subroutine  any  number  of  times,  and  we  show  some 
of  the  results  for  up  to  3  cycles.  The  best  results  were  obtained  by 
using  two  cycles  with  k  =  1.  However,  this  equalization  of  ND(1)  and 
ND(2)  is  purchased  at  the  cost  of  a  considerable  reduction  in  the  number 
of  data  points,  which  adversely  affects  the  performance  of  the  algorithms. 
The  data  set  F  was  obtained  by  numerically  raising  the  threshold  level  of 
Radar  2.  It  turned  out  that  the  disparity  was  almost  completely  reversed 
at  the  lowest  possible  threshold  setting,  so  that  the  goal  of  equalizing 
the  number  of  detections  could  not  be  attained  in  this  way.  In  addition, 
we  also  examined  the  results  of  adjusting  the  threshold  levels  of  both 
radars;  however,  we  found  that  the  accuracy  of  the  corresponding  crab 
angle  estimates  was  not  improved,  and  these  results  are  not  shown.  In  the 
sequel  we  shall  show  the  results  of  applying  various  cross  correlation 
techniques  to  the  two  sets  A  and  F,  and  of  using  various  centroid  type 
estimates  on  the  six  sets  A-F. 


TABLE  2 

Numbers  of  Detections 

ABC 

Original  Data  Outliers  Removed  Outliers  Removed 

L  =  y  +  a;  1  cycle  L  *  p  +  2o;  1  cycle 


ND(1) 

ND(  2) 

RATIO 

ND(  1) 

ND(2) 

RATIO 

ND(  1) 

ND(2) 

RATIO 

File 

3 

747 

855 

1.14 

497 

517 

1.04 

723 

839 

1.16 

File 

5 

665 

876 

1.32 

429 

549 

1.28 

641 

859 

1.34 

File 

6 

776 

889 

1.15 

523 

557 

1.07 

741 

868 

1.17 

File 

7 

827 

921 

1.11 

544 

571 

1.05 

799 

908 

1.14 

D 

E 

F 

Outliers  Removed 

Outliers  Removed 

Threshold  Adjust 

L  =  y 

+  a;  2  cycles 

L  -  y 

+  a;  3  cycles 

ND(  1) 

ND(2) 

RATIO 

ND(1) 

ND(2) 

RATIO 

ND(  1) 

ND(  2) 

RATIO 

File 

3 

295 

294 

1.00 

169 

160 

0.95 

747 

626 

0.84 

File 

5 

276 

312 

1.13 

165 

184 

1.12 

665 

681 

1.02 

File 

6 

313 

324 

1.04 

183 

180 

0.98 

776 

667 

0.86 

File 

7 

324 

328 

1.01 

179 

186 

1.04 

827 

687 

0.83 
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Another  indicator  of  the  dissimilarity  between  the  two  records  is  pro¬ 
vided  by  the  dissimilarity  between  their  autocorrelation  functions,  each  of 
which  is  obtained  by  cross  correlating  a  record  with  a  copy  of  itself  at 
lagged  values  of  N  and  R.  The  autocorrelation  functions  of  AMP^(N,R)  and 
AMP2(N,R)  are  shown  in  Figures  1  and  2.  The  data  is  from  File  5.  The 
entries  are  the  values  of  the  correlation  coefficients  multiplied  by  a 
factor  of  100,  except  for  the  first  column  which  shows  the  lagged  value  AN 
of  N.  The  lagged  values  of  AR  increases  as  one  moves  horizontally  from 
left  to  right,  but  their  values  are  not  shown  since  they  are  not  involved 
in  the  calculations  of  the  At  used  in  Eq.  (1).  The  formulas  used  to  calcu¬ 
late  the  correlation  functions  and  centroids  are  given  in  the  Appendix. 

The  peak  value  of  an  autocorrelation  is  unity  (by  construction),  and  the 
next  largest  values  are  p  =  .42  in  Figure  1  and  p  3  .20  in  Figure  2.  This 
difference  is  significant. 

The  cross  correlation  function  of  AMPj  with  AMP2  is  shown  in  Figure 

3.  The  peak  value  p  =  .14  occurs  at  AN  =  50,  and  is  marked  with  an 
asterisk.  The  small  size  of  this  peak  value  is  another  indicator  of  the 
dissimilarity  between  the  two  records. 

The  column  containing  the  peak  value  of  p  in  Figure  3  displays  a  lobed 
structure  which  is  very  flat  and  noisy  in  appearance,  which  makes  it 
unsuitable  for  accurate  peak  location.  For  this  reason  we  applied  various 
smoothing  operations  to  the  data,  which  are  described  in  the  next  section. 

4.  RESULTS  USING  CROSS  CORRELATION 

Two  types  of  data  smoothing  schemes  were  used  in  applying  the  cross 
correlation  technique.  First,  the  (N,R)-space  was  divided  into  rectangular 
blocks  consisting  of  n  successive  pulses  and  r  contiguous  range  bins,  and 
the  amplitudes  in  each  block  were  summed  to  form  two  new  arrays  which  are 
then  cross  correlated.  The  advantage  of  this  technique  is  that  it  produces 
a  correlation  function  whose  lobes  are  sharper  and  smoother  than  those 
obtained  from  the  original  data.  The  disadvantage  is  that  it  introduces  a 
quantization  error  which  increases  with  increasing  n.  Second,  the  corre¬ 
lation  function  itself  was  smoothed  by  applying  a  5  point  moving  average 
along  each  column,  (i.e.,  fixed  range  lag). 

The  results  for  n  =  5  and  r  *  3  are  shown  in  Table  3.  The  aircraft 
velocity  was  approximately  250  ft. /sec,  and  referring  to  Eq.  (1)  with  the 
indicated  value  of  D,  it  follows  that  the  quantization  error  for  n  =  5  is 
approximately  0.34°.  Other  values  of  n  and  r  sometimes  produced  a  better 
lobed  structure,  but  the  corresponding  angle  estimates  were  never  satis¬ 
factory.  Three  different  values  are  tabulated,  corresponding  to  three 
different  methods  for  peak  location: 

(M).  The  location  at  which  p  is  maximum. 

(MAM).  The  location  at  which  the  5  point  moving  average  is  maximum. 

(BS).  The  result  obtained  by  beam  splitting  the  5  point  moving 
average  by  taking  the  average  of  the  two  3  dB  points. 

For  convenience,  the  "true"  values  of  apparent  crab  angle  are  shown  in 
parentheses;  i.e.,  from  the  "sum"  column  of  Table  2.  Table  4  shows  the 
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results  when  the  same  technique  is  applied  to  data  set  obtained  by  applying 
the  threshold  adjustment  discussed  above  in  constructing  the  data  set  "F" 
of  Table  2.  The  results  are  exactly  the  same  except  at  two  places. 


TABLE  3  TABLE  4 


Cross  Correlation  Estimates  Estimates  with  Threshold  Adjust 


M 

MAM 

BS 

M 

MAM 

BS 

File 

3 

(-0.25) 

0° 

+0.37 

+0.56° 

File 

3 

0° 

+0.37° 

+0.56 

File 

5 

(+5.3) 

-1.48 

-1.48 

+2.22 

File 

5 

-1.48 

-1.48 

+2.22 

File 

6 

(+1.1) 

+0.68 

+1.03 

+1.03 

File 

6 

+0.68 

+1.03 

+1.03 

File 

7 

(-1.4) 

+0.35 

+0.35 

-0.35 

File 

7 

+0.35 

0 

-0.35 

File 

4 

(-2.6) 

-0.10 

-0.64 

-0.96 

File 

4 

-0.10 

-0.64 

-0.96 

File 

9 

(-6.0) 

+8.17 

+8.17 

+8.46 

File 

9 

+8.17 

+8.17 

+8.46 

File 

10  (+2.3) 

+2.28 

+2.59 

+1.66 

File 

10 

+2.28 

+2.59 

+1.97 

Among  the  first  four  "initial"  files,  the  results  are  always  bad  for 
the  large  crab  in  File  5,  but  marginally  acceptable  results  are  contained 
in  the  BS  column  for  Files  3,  6,  and  ?.  The  optical  data  from  File  5  was 
examined  very  closely,  and  the  bad  results  obtained  from  the  "additional" 
Files  4  and  9  show  that  File  5  was  not  an  anomaly.  The  result  for  File  9 
is  particularly  bad  because  of  the  sign  reversal. 

Tables  5  and  6  show  the  results  of  applying  the  same  calculations  to 
"binary  weighted"  data,  i.e.,  data  for  which  AMPj(N,R)  is  replaced  by  0 
or  1  depending  on  whether  a  detection  occurs  at  the  indicated  value  of  N 
and  R.  This  technique  produced  an  increase  in  the  peak  values  of  p,  but 
did  not  improve  the  crab  angle  estimates.  The  peak  values  of  p  are  shown 
in  Table  7,  for  both  the  ordinary  and  binary  weighted  data.  The  simple 
maxima  of  the  cross  correlation  functions  are  denoted  by  "M",  and  the 
maxima  of  the  5  point  moving  average  by  "MAM".  The  effects  of  threshold 
adjustment  on  the  peak  values  of  p  were  miniscule,  and  are  not  shown. 

The  peak  values  of  p  for  Files  4,  9,  and  10  (with  "ordinary"  data)  are 
.44,  .43,  and  .40,  respectively,  and  the  observed  tendency  for  the  peak 
values  of  p  to  decrease  as  the  magnitude  of  the  crab  angle  increases  is 
very  slight,  probably  due  to  the  complicating  effects  of  range  differ¬ 
entials  and  aircraft  roll. 


TABLE  5 


TABLE  6 


Cross  Correlation  Estimates  Cross  Correlation  Estimates 

of  Crab  Angle  with  of  Crab  Angle  Data  with  Binary 

Binary  Data  Data  and  Threshold  Adjust 


M 

MAM 

BS 

M 

MAM 

BS 

File 

3 

(-0.25) 

+0.56° 

+0.56° 

+1.85° 

File 

3 

+0.74° 

+0.74° 

+1.48° 

File 

5 

(+5.3) 

-1.85 

-1.85 

-0.37 

File 

5 

-1.85 

-1.85 

-0.56 

File 

6 

(+1.1) 

+0.85 

+0.34 

-0.51 

File 

6 

+0.85 

+0.34 

-0.34 

File 

7 

(-1.4) 

+0.35 

+0.69 

-0.17 

File 

7 

+0.35 

+0.35 

-0.17 

TABLE  7 

Peak  Values  of  Correlation  Coefficient  p 

ORDINARY  DATA  BINARY  DATA 


M 

MAM 

M 

MAM 

File 

3 

.44 

.41 

.51 

.50 

File 

5 

.36 

.33 

.44 

.42 

File 

6 

.48 

.45 

.57 

.56 

File 

7 

.52 

.49 

.60 

.59 

Figure  4  shows  the  autocorrelation  function  used  in  obtaining  the 
results  for  File  5  in  Table  3.  Asterisks  mark  the  locations  of  the  peak 
and  the  two  3  dB  points  of  the  5  point  moving  average.  The  maximum  value 
in  an  adjoining  column  was  very  close  to  the  peak,  and  is  marked  with  a 
solid  disk.  Its  location  is  closer  to  that  corresponding  to  the  actual 
crab  angle.  By  reason  of  the  quantizing  of  the  (N,R)-space  into  5x3 
rectangular  blocks,  each  unit  increment  of  lag  AN  corresponds  to  a  time 
delay  At  of  5  ms.  Therefore  the  "true"  peak  location  occurs  at  AN  =  15. 
Figure  5  shows  the  one-dimensional  slice  of  the  correlation  function 
taken  along  the  main  column  (containing  the  peak).  The  estimates  for 
File  5  were  very  bad  and  Figure  5  shows  how  this  happened:  The  correla¬ 
tion  function  for  this  particular  file  is  non-symmetric  and  has  multiple 
peaks. 

Numerous  other  cross  correlation  estimators  were  examined  besides 
those  described  above,  viz.,  (1)  different  values  for  the  parameters 
(n,r)  in  the  (N,R)-space  quantization,  (2)  different  thresholds  applied 
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to  one  or  both  radars,  (3)  logarithmic  weighting  (replacing  AMPi(N,R)  with 
its  logarithm),  (4)  one  dimensional  correlation  obtained  by  weighting  each 
pulse  by  the  number  of  detections  or  sum  of  all  values  (setting  r  =  ®  in 
(1)).  These  four  techniques  can  be  combined  in  numerous  ways,  but  none  of 
the  dozens  examined  produced  acceptable  results. 

5.  CENTROIDS 


The  centroids  (or  means)  and  112  o£  the  arrays  AMPj  and  AMP2  were 
calculated  according  to  formulas  given  in  the  appendix,  and  the  correspond¬ 
ing  crab  angle  estimates  were  the  obtained  by  setting  At  =  (^2  ~  yj.)/1000 
in  Eq.  (1).  Three  types  of  weighting  scheme  were  used  in  the  calculation 
of  centroids: 

Type  1.  Each  pulse  weighted  by  the  number  of  detections 

Type  2.  Each  pulse  weighted  by  0  or  1  depending  on  whether  a 
detection  occurred  in  that  pulse. 

Type  3.  Each  detection  weighted  by  AMPj(N,R). 

Type  3  is  the  more  conventionally  defined  centroid.  In  addition  to  these 
three  weighting  schemes,  there  are  six  categories  A-F  which  correspond  to 
those  in  Table  2.  We  recall  that  A  is  obtained  from  the  original  untouched 
data,  B-E  from  data  with  outliers  removed,  and  F  from  data  which  has  been 
threshold  adjusted  as  discussed  in  Section  3.  There  are  therefore  18 
different  kinds  of  centroids.  The  results  are  shown  in  Table  8.  The 
results  for  the  "additional"  Files  4,  9,  and  10  were  only  calculated  for 
the  categories  A,  B,  C. 

Again,  other  types  of  centroids  were  calculated  by  using  logarithmic 
weights,  by  adjusting  the  thresholds  on  one  or  both  radars,  and  by  varying 
the  parameters  in  the  outlier  removal  scheme.  But,  as  before,  acceptable 
results  were  never  achieved. 

6.  SUMMARY  AND  CONCLUSIONS 


The  system  concept  is  based  on  the  idea  that  the  returns  from  the  two 
radars  should  be  similar  except  for  the  time  lag  produced  by  the  crab; 
however,  the  two  radar  returns  were  found  to  be  dissimilar.  One  radar 
always  had  more  detections  than  the  other,  the  two  autocorrelation  func¬ 
tions  were  different,  and  the  peak  values  of  the  cross  correlation  were 
small.  This  dissimilarity  between  the  two  records  was  probably  caused  by 
range  differentials,  aircraft  roll,  and  the  crab  angle  itself. 

In  addition  to  the  examination  of  cross  correlation  techniques,  we  also 
considered  the  use  of  estimates  derived  from  the  centroids,  which  are  much 
easier  to  implement  numerically. 

Hundreds  of  different  types  of  crab  angle  estimators  were  examined, 
each  one  corresponding  to  different  possibilities  for  thresholding,  data 
weighting,  data  smoothing,  and  outlier  removal.  None  of  these  techniques 
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Estimates  Derived  from  Centroids 
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produced  acceptable  results,  and  the  performance  at  the  two  largest  crab 
angles  (5°  and  6°)  was  always  extremely  poor.  We  therefore  conclude 
that  any  further  examination  of  this  system  would  not  be  warranted. 
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Figure  4.  Cross  Correlation  Function  for  the  Quantized  Data  from  Fixe  5 

n  =  5  r  =  3 
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APPENDIX 


The  cross  correlation  coefficient  p(AN,  AR)  was  computed  according  to 
N2  R2  r  -i 


p(AN, AR)= 


m-  -» 

Y  Y  AMP  (N,R)-AMP  (N+AN,  R+AR) 
N=N1  R=R1  L  1  Z  J 


'•  (A-l) 


N2 

R2  r  19  1'2 

N2 

R2  r  1 2 

E 

E  AMP.(N,R) 

E 

Y  AMP  (N+AN,  R+AR) 

N=N1 

R=R1  L  J 

N=N1 

R=R1  L  J 

The  limits  of  summation  Nl,  N2,  R1  and  R2  vary  with  AN  and  AR,  and  are 
chosen  so  that,  in  effect,  the  data  is  zeroed  where  the  two  records  do  not 
overlap.  In  the  calculation  of  p  with  binary  weighted  data,  we  merely  set 
AMPi  and  AMP2  equal  to  0  or  1  depending  on  whether  a  detection  occurs  at 
the  indicated  value  of  the  arguments.  It  is  a  common  practice  to  subtract 
the  means  from  the  amplitudes  in  this  calculation.  However,  this  procedure 
does  not  affect  the  location  of  the  peak. 

The  two  centroids  and  P2  of  Type  1  were  calculated  according  to 

E  N  •  V  (N) 

P1  =  ^ -  *  (A"2> 

E  vjoo 


where  vj(N)  is  the  number  of  detections  at  the  N'th  pulse  of  the  I’th  radar. 
The  centroids  of  Type  2  are  given  by 


E  N  •  ♦  (N) 

=  _N _ 

1  E  <MN> 

11  ^ 


(A- 3) 


where  <f>j(N)  is  0  or  1  depending  on  whether  a  detection  occurs  at  the  N'th 
pulse.  The  more  conventionally  defined  centroids  of  Type  3  are  given  by 


EE  N  •  AMP  (N,R) 
N  R 

j  =  - - — ■ — — 

1  EE  AMPjCN.R) 


(A- A) 


Note  that  the  Type  1  centroids  are  equivalent  to  Type  3  centroids  with 
binary  weighted  data. 
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